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Abstract. Data on stopping in intermediate-energy central heavy-ion collisions are analyzed following transport theory based 
on the Boltzmann equation. In consequence, values of nuclear shear viscosity are inferred. The inferred values are significantly 
larger than obtained for free nucleon dispersion relations and free nucleon-nucleon cross sections. 
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INTRODUCTION 

In central reactions of heavy nuclei, momentum is transferred between matter originating from opposing nuclei. 
Towards the end of a reaction, the matter can be described as locally equilibrated, in terms of local temperature field 
folded with a field of local collective velocity. Obviously, in a reaction, dissipation takes place. Looking at appropriate 
data from the reactions, one can ask about the pace of dissipation and examine what that pace tells about the general 
dissipative properties of nuclear matter. 

A theoretical model for the central reactions needs to be capable of describing different stages of a reaction, 
nonequilibrium and equilibrium. Here, we shall rely on a single-particle description, in terms of the nucleon Wigner 
functions f(p,r,t) obeying a set of Boltzmann equations: 
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In the above, the single-particle energy £ is a variational derivative of the energy represented as functional of 
the Wigner functions, e = (27r)~ 3 dE/df. The l.h.s. of (Q~|) accounts for changes in the Wigner function due to the 
movement of particles, at velocity v = de/dp, and due to their acceleration on account of the single-particle energy 
changing with position, with F = —de/dr representing a force on the particle. The r.h.s. of the equation accounts 
for changes in / due to collisions. The pace of approach to local equilibrium is governed by cross-sections a, most 
often assumed to coincide with those in free space. The factors / are Pauli-blocking factors, / = 1 — /. In models of 
reactions where the Boltzmann set is not directly followed, it is still common to incorporate elements of the Boltzmann 
equation, such as the scattering governed by cross sections J2]. 

In comparing a transport model to reaction data, assumptions within the model are adjusted until a reasonable 
agreement reached. Obviously, some data test some assumptions better and some worse. An issue is the universality 
of drawn conclusions. Information pertaining to a transient reaction stage alone or, even worse, just to specific model, 
can be of very limited utility. For that reason, while a model is used to describe a reaction through its nonequilibrium 
stages, the same model is typically used to extrapolate the conclusions to those pertaining to finite- or zero-temperature 
equilibrium, and referring then to energy, pressure or optical potential JH S S S lal • Here, see also [7, 8], we shall 
consider characterization of the nuclear system for weak deviations from equilibrium, in terms of shear viscosity 
coefficient. Just as the equation of state, the macroscopic transport coefficients might be assigned to nuclear matter even 
for models partly or fully phenomenological in nature, provided one could legitimately tie the variation in transport 
properties to particular observables. 

The Boltzmann set (Q~|) should apply to nuclear systems at low nucleon density n combined with moderate to high 
temperatures T. At high densities, when the binary-collision regions would likely overlap, the equation provides 
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FIGURE 1. Schematic illustration for momentum transport in a reaction, between the projectile and target zones. In a semipe- 
ripheral reaction, the collective velocity has predominantly longitudinal z-components and changes in transverse directions such 
as the x-direction along the reaction plane. Associated with the velocity gradient, close to equilibrium, is the transport of the z- 
component of momentum in the x direction. Any location in the matter is reached by nucleons from about one mean-free-path X 
away, contributing to the momentum flux. 



just phenomenological description. However, when the collisions are frequent in a system, the system approaches 
a local thermal equilibrium and begins to behave hydrodynamically. In the latter limit, microscopic details behind the 
hydrodynamic behavior might not be important, as long as macroscopic properties are properly reproduced. 



For weak gradients in a largely equilibrated system, fluxes of macroscopic quantities, leading to dissipation, are 
proportional to gradients within the system. The Curie principle, stating that transformation properties for interrelated 
fluxes and gradients must be the same, allows for sorting out the possible relations. The shear viscosity coefficient tj, 
that we shall be after, is the coefficient of proportionality between anisotropy of momentum-flux tensor, inducing 
dissipation, and velocity gradients. When local velocity u is directed along a specific z-direction, as is the case 
approximately in a semiperipheral nuclear reaction, see Fig. [T] and when u changes in value in the perpendicular 
direction x, the flux IF-*, of z-component of momentum in x-direction, is 



to linear order in gradients. 

A simple estimate for the viscosity coefficient may be arrived at by carrying out mean-free-path considerations. 
Thus, for transport in the x-direction, different locations along the x-axis will matter. Any position along the x-axis, 
see Fig.Q] will be reached by particles that start out one mean free path X away. Only about 1/6 of particles at the 
starting location move, at typical speed v^n, towards the point of interest. The average z-momentum, they bring in, is 
mu z , where u z is that for their starting location. Accounting for the fluxes due to particles moving up and down the 
x-axis yields 
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where n is particle density. This yields then an estimate for the viscosity coefficient of 
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where, in the quantitative estimate, we have used values representative for reactions at moderate beam energies. 
Proportionality of the viscosity to the mean free path implies that the viscosity is inversely proportional to the 
interaction cross sections, or otherwise exhibits a negative correlation with stopping, since X ~ 1 /{no). 



When comparing data to the results of the Boltzmann set, to infer nuclear viscosity, one needs to determine first 
the viscosity for the set. For this, one must consider a system governed by £[), close to equilibrium. The r.h.s. of the 
equations in the set vanishes for local equilibrium distributions, 

f e i = 1 (5) 
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where u, /X, and T are equilibrium parameters that depend on position. However, in the latter case, the distributions 
cannot solve the set, because of the derivatives evaluated on the l.h.s. With (Q}, the local equilibrium distributions (0 
need to be to corrected to at least the first order in gradients in equilibrium parameters, for these distributions to solve 
the Boltzamnn set. Notably, the local-equilibrium distributions <(5j would yield vanishing dissipative contributions to 
macroscopic fluxes, in particular to the anisotropy of the tensor of momentum flux. On the other hand, corrections 
linear in gradients will yield contributions to fluxes that are linear in gradients, that are of interest in the context 
of ©. Following the Chapman-Enskog method, one can, in fact, systematically seek a solution to ([1]) around local 
equilibrium, by expanding the distributions in derivatives 

f i =tf ) +fV+f? ) + ... (6) 

where fj- ' = f eq and is of the n'th order in gradients. The n'th order corrections may be obtained from the 
Boltzmann set by inserting the (n — l)'th terms into the l.h.s. 

For inferring shear viscosity, the terms are important. Upon inserting ff q to the l.h.s. of (fl}, the form of fj- 
emerges lr9l[Tofl: 

tf^tifPil-fh (7) 
where 0, are smooth functions of position and momenta. At low temperatures, the distributions are modified, com- 
pared to local equilibrium, in the region of Fermi surface. Following the Curie principle, anisotropy of symmetric 
momentum-flux tensor should be driven by the anisotropy of symmetric tensor of velocity gradient: 

0; = b t [pkPt - y §ki^J fvjt u e + W e u k - j 5u V„u„^j , (8) 

where bt are, generally, functions of momentum magnitude. Assuming that b, change weakly momentum, one can 
arrive at a closed expression J*! [T/J] for the viscosity, accurate in practice down to few percent for a given cross 
section, and provided here in the simplified form for symmetric matter: 

n = s S a^ 2 ) 2 (9) 
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The shear viscosity coefficient indeed comes out in (0 inversely proportional to the cross sections. However, the 
cross-sections in the viscosity get weighted with Pauli-blocking factors. Additional weighting there indicates that 
the cross sections at high relative-momenta and large scattering-angles matter more for the viscosity than the cross 
sections at low momenta and small angles. The left panel in Fig.|2]shows the shear viscosity from Boltzmann equation, 
calculated with cross sections and velocities such as in free space. For typical conditions in a reaction, at moderate 
incident energies, the calculated viscosities are higher than in the simple mean-free-path estimate. At low temperatures, 
the viscosity values diverge due to Pauli blocking of collisions. 

We shall next examine whether data pertinent to dissipation of momentum lfl"ll[l2tl justify the use of free elementary 
cross sections in the description of nuclear dynamics. Other types of data on nuclear dissipation, of general interest 
especially in the context of exotic beams, are those on equilibration of neutron-proton asymmetry ]9l 0, [T4IJT5 
Previous data analyses have established a lowering of nucleon effective masses with increase in nuclear density |l| 

DATA COMPARISONS 
vartl Observable 



The first data, that we are going to compare our theory to, are those on the degree of isotropization of momentum 
distributions in symmetric central collisions. To quantify the degree of isotropy in the final-state of a reaction, the 
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FIGURE 2. Viscosity in symmetric nuclear matter at different densities, as a function of temperature. Left panel shows viscosity 
calculated with free-space cross-sections and velocities. Right panel shows viscosity calculated with in-medium cross-sections and 
effective masses adjusted to heavy-ion data. 



FOPI Collaboration [11] has introduced the observable 

Ay, 

vartl = — , (10) 
Ay/ 

that is the ratio of the widths of rapidity distributions in the transverse and longitudinal directions. The rapidity for the 
longitudinal direction is defined in the standard manner. The rapidity for the transverse direction is defined by replacing 
the beam direction, in the standard definition, by a random transverse direction. Reaching isotropy in a reaction would 
produce vartl ~ 1. Transparency effects would yield vartl < 1. Finally, a strong hydrodynamic behavior might yield 
vartl > 1 in a central reaction. As systems should evolve, from transparency towards a hydrodynamics behavior, with 
increasing system mass, vartl is expected to grow with system mass. 

Excitation functions for central Au + Au and Ca + Ca collisions, from measurements of the FOPI Collaboration Hill , 
are shown in Fig. [3] The FOPI Collaboration has calculated the vartl values using all fragments with Z < 10. To 
illustrate the effect of including the fragments with relatively high-Z values, we include also an Au + Au result 
from the INDRA Collaboration 11611 . obtained including fragments with Z = 1 and 2 only. In central collisions, the 
intermediate-mass-fragments with Z > 3 are abundant at beam energies around 100 MeV/nucleon, but become less 
frequent at 400 MeV/nucleon and above. Also, such fragments are more frequent in a heavy system such as Au + Au 
than the relatively light Ca + Ca. 

As is apparent in Fig. [3] when looking at the Z < 10 fragments, the proximity to isotropy is observed only at 
intermediate energies in Au + Au collisions, with vartl > 0.85 there at 400 — 800 MeV/nucleon. At low energies, the 
original nuclei might not interpenetrate enough for the momentum dissipation to get completed and Pauli principle 
might further suppress collisions. At high energies, elementary cross sections become forward peaked, leading to 
transparency. 

To facilitate comparisons to data, we include in the model calculations, at low densities, the production of deuterons, 
tritons and helions in few-nucleon collisions lfl7tl . Since, however, the model does not predict yields of heavier 
fragments, the comparisons to the FOPI data will not be meaningful below 250 MeV/nucleon for Au + Au reactions. 
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FIGURE 3. Excitation function of the vartl observable in central Au + Au (left panel) and Ca + Ca (right panel) collisions. 
Symbols represent data of the FOPI fTHl (triangles) and INDRA 1 16] (square) Collaborations. The FOPI results are obtained using 
all fragments with charge number Z < 10. The INDRA result is for Z = 1 and 2. Lines represent results from the Boltzmann 
equation model including fragments with mass number A < 3, for different assumptions on elastic elementary cross sections. 



Also, in comparisons, we will manipulate only elastic elementary cross sections only, dominating reactions below 
800 MeV/nucleons. Thus, most relevant conclusions will pertain to the energy region of 250 — 800MeV/nucleon. 

It is apparent in the left panel of Fig. [3] that the model calculations with free cross sections grossly overestimate 
stopping data, yielding even values of vartl > 1 below 800 MeV/nucleon, in the indication of a strong hydrodynamic 
behavior. Clearly, the free cross sections produce excessive stopping in those reactions. Subsequent issue is of how the 
cross sections should be modified and, more accurately, reduced in the medium. There had been calculations in the 
literature, carried out for equilibrated nuclear-matter, of changes in the cross sections due to effects of effective mass 
and effects of the Pauli principle on intermediate states in two-nucleon scattering lfl8l [l9l l20ll . General finding has 
been that of lowering of the low-energy nucleon-nucleon cross sections in the medium. Testing those cross sections, 
we carry out Boltzmann-equation simulations employing parameterized results on the in-medium reduction of the 
cross-sections, by the Rostock group and their collaborators lfl8l[T9tl . The vartl results shown in Fig. [3] are reduced 
compared to those for the free cross-sections, but they are still excessive compared to the data. 

The nuclear-matter calculations of cross sections, such as 11181,1 19112011 . do not account for overlapping of the collision 
regions at high densities, that compete against each other. To account phenomenologically for a unitary saturation 
taking place, we require that the nucleon-nucleon cross section cannot exceed a size imposed by interparticle distances 

0-<0-o = vn~ 2/3 , (11) 
where v is a parameter of the order of 1 . To realize this limit in practice, we parameterize the in-medium cross with 
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where v is adjusted. In the low-density limit, n — > 0, the in-medium cross-section approaches then the free cross- 
section, o — > Of re e- At high-density, n — > oo, the in-medium cross section approaches Oo from below, o /* Oo- From the 
Au + Au results for different v in Fig. [3] it is apparent that v = 0.6 is closest to the data in the relevant energy region. 

If we next turn to the Ca + Ca vartl results in Fig. [3] we find again that calculations with free cross sections 
again strongly overestimate the measurements. However, the v = 0.6 calculations, favored in the Au + Au case, also 
overstimate the data. The Ca + Ca data rather favor v = 0.4. However, an issue in the Ca + Ca reactions is of a less 
precise determination of the centrality of reactions than in Au + Au. Lower multiplicities in Ca + Ca collisions, than 
in Au + Au, produce fluctuations for observables used to constrain reaction centrality. If we assume that the collisions 
in the measurements have not been fully central, we can get agreement between calculations and data for V ~ 0.6. 




FIGURE 4. Schematic illustration of measurements of linear momentum transfer. Among products of the reaction, the most 
massive fragment is identified. Its longitudinal velocity is compared to the center of mass velocity for the system. 
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FIGURE 5. Excitation energy for the measure of stopping (vii/v C m) in central Ar + Cu (left panel) and Ar + Ag (right panel) 
collisions. The velocity component vm is that of the heaviest fragment emitted from the reaction. Symbols represent data of Ref. ITHl 
and lines represent results of Boltzmann-equation simulations employing different assumptions on elementary in-medium cross 
sections. 

Linear Momentum Transfer 

Another observable, pertinent to stopping, has been the linear momentum transfer (LMT), used, in particular, for 
quantifying central asymmetric collisions of heavy ions. In the measurements [12], the heaviest emitted fragment has 
been identified, see Fig. [4] Most likely that fragment originates from the target residue left after the fast initial stage 
of the reaction. The average longitudinal component of the fragment velocity should reflect the average velocity of 
the residue. A scale-invariant measure of the degree of stopping in a reaction may be obtained by taking a ratio of the 
average fragment velocity to the center-of-mass velocity for the system as a whole. If the target and projectile fuse, 
then the ratio should be close to 1, (vi|}/v C m — 1. Little stopping should be characterized by (vi|)/v cm <C 1. 

Figure[5]shows measured and calculated excitation functions for (vii ) / v em in central Ar + Ag and Ar + Cu collisions. 
At low beam energies, the nuclei appear to fuse. As beam energy increases stronger and stronger transparency sets 
in. Calculations with free cross definitely overestimate the stopping in the reactions. The Rostock in-medium cross 
sections reduce stopping, but not quite enough. For the cross-sections of Eqs. (fT2l i and (fTTT i. getting close to the data 
requires going down in v to the vicinity of v = 0.6. 

NATURE OF CONCLUSIONS 

It is apparent that the use of free cross sections in the Boltzmann-equation simulations yields too much stopping. 
The nature of the implied in-medium reduction of cross-sections, however, is not obvious. Specifically, different types 
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FIGURE 6. Measures of elementary collisions in the simulations of central Ar + Ag reaction at 90 MeV/nucleon, for different 
assumptions on elementary cross sections, as a function of time t in the reaction. The left panel shows net collision number. The right 
panel shows the number of collisions weighted with the weight p\ 2 sin 8 2 that is applied to the collision rate in the expression l|9j 
for viscosity. 



of cross-section reduction can produce similar stopping. Thus, with v ~ 1 in Eqs. ( fTTT ) and ( fT2b . we can get nearly 
identical stopping results in terms of vartl and (v||/v cm ), for the reactions in Figs.|3]and|5] However, collision counts 
for the simulations utilizing the two cross-section reductions are then vastly different, see the left panel in Fig.|6]for 
results for the Ar + Ag reaction at 90 MeV/nucleon. Compared to the simulation with free cross sections, the collision 
number drops just by 25% and by a factor of 4, respectively, for the Rostock and v = 1 cross sections, respectively! 

As we have stated, however, it is not clear to what extent the separate nucleon-nucleon collisions can be, in reality, 
identified within a reaction. Even if they can be identified, some, especially those at low relative momenta and those 
with forward scattering angles, may matter little for the reaction dynamics. Around equilibrium, solely macroscopic 
properties of matter would be important for the dynamics and the same dynamics would result from different cross 
sections as long as those cross sections gave the same transport coefficients. For dissipation of momentum, the 
dominant role would be played by the viscosity coefficient and, in fact, that coefficient might not even be tied to 
a medium for which the kinetic limit applies. In the coefficient for the kinetic limit, Eq. (O, the cross section is 
multiplied by the weight p\ 7 sin 9 2 , which emphasizes collisions at high relative momentum, populating wide angles. 
Collisions weighted in this fashion should matter for dissipation of momentum close to equilibrium. The right panel 
in Fig.|6]next shows the collision numbers for different cross sections weighted with the viscous weight. The weighted 
collision numbers are now similar for the two in-medium cross-sections that produce similar stopping, consistently 
with the discussion above. With a similar count of the weighted collisions, the two cross would produce similar shear 
viscosity coefficients for conditions such as in the reactions. 

The right panel in Fig.[2]shows the shear viscosity coefficient calculated with the v = 0.6 cross-section values, that 
yield stopping of such order as observed, and with in-medium dispersion relations established earlier [1]. Both the 
lowering of cross sections and of effective masses contribute to an enhancement of the coefficient values, compared to 
the results obtained when disregarding in-medium effects, in the left panel of Fig. [2] 



SUMMARY 

We have confronted the measurements of stopping in central nuclear collisions, in terms of the vartl observable and 
the linear-momentum transfer, with the predictions of stopping from a Boltzmann-equation set. Predictions utilizing 
free nucleon-nucleon cross-sections strongly overestimate the stopping, irrespectively of the considered observable or 
the system. Predictions utilizing cross-sections from nuclear-matter calculations 11181 [l9l 12011 yield less stopping than 
for the free cross sections, but still too much compared to data. The stopping data alone do not constrain details in 
the in-medium cross-section alone. Similar stopping predictions can be arrived at within calculations where collision 
numbers differ by a factor of 3. Fortunately, for drawing conclusions from the stopping data, similar stopping results 



are correlated to similar predictions for nuclear shear viscosity. With inclusion of the in-medium effects, the deduced 
viscosity values are significantly larger than anticipated in the absence of such effects. Systematics of the nuclear 
viscosity values is, in particular, of interest in assessing how a nuclear system may approach the limit of near-perfect 
liquid with increase of temperature [21]. The consideration of the latter limit involves also analysis of entropy values 
that may be deduced from fragment yields in the reactions. 
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